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Abstract 


Sodium borohydride (NaBH) hydrogen storage systems offer many advantages for hydrogen storage applications. The physical processes 
inside a NaBH; packed bed reactor involve multi-component and multi-phase flow and multi-mode heat and mass transfer. These processes are 
also coupled with reaction kinetics. To guide reactor design and optimization, a reactor model involving all of these processes is desired. A one- 
dimensional numerical model in conjunction with the assumption of homogeneous catalysis is developed in this study. Two submodels have been 
created to simulate non-isothermal water evaporation processes and pressure drop of two-phase flow through the porous medium. The diffusion 
coefficient of liquid inside the porous catalyst pellets and the mass transfer coefficient of water vapor are estimated by fitting experimental data 
at one specified condition and have been verified at other conditions. The predicted temperature profiles, fuel conversion, relative humidity and 


pressure drops match experimental data reasonably well. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


A sodium borohydride hydrogen generator is unique because 
both reactants can be stored together, and hydrogen is generated 
by passing sodium borohydride solution through a catalyst bed 
to initiate hydrolysis reaction as [1]: 


NaBH; + 2H,0*“S"NaBO, + 4H> + Heat (1) 


The effects of catalysts, pH and temperature on sodium boro- 
hydride hydrolysis reaction were discussed in Ref. [2]. Our 
experimental paper [3] discussed system-level experiments on a 
1 kW, sodium borohydride hydrogen generator and exposes sol- 
ubility issues that may limit the maximum usable concentration 
to approximately 15%, which may preclude automotive appli- 
cations; nevertheless, sodium borohydride systems may still 
find applications in portable electronic devices and other niche 
areas [4,5]. Most prior work on sodium borohydride systems has 
focused on experimental testing, and no work on system-level 
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reactor modeling has been reported to date. The processes in 
the reactor are quite complex, involving multiple components 
(NaBH4, NaOH, NaBO2, H20, H2) and multiple phases (liq- 
uid and gas). In addition, general liquid phase reactions inside a 
packed bed reactor accompanied by significant water evapora- 
tion have not received attention in the literature. As a result, 
a significant need exists for a sodium borohydride hydroly- 
sis reactor model to enable reactor design and optimization. 
The development of such a reactor model will also facilitate 
the study of hydrogen storage systems using other chemical 
hydrides. Thus motivated, we have developed a one-dimensional 
numerical model in conjunction with the assumption of homoge- 
neous catalysis and have validated this model with experimental 
data. 


2. Experiments 


A 1kW, sodium borohydride hydrogen generation system 
(by assuming a fuel cell efficiency of 50%) has been estab- 
lished for system-level studies. The 1 kW. hydrogen generation 
apparatus was described in detail in a previous paper [3] and 
is only summarized here. Fig. 1 shows the section view of the 
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Nomenclature 

ag external surface area per volume of catalytic bed 
(m? m~?) =6(1 — £)/dp for packed bed (m? m~*) 

At cross-sectional area of the reactor (m?) 

Cp,f specific heat of the fuel (kJ kg~! k~!) 

Ca molar concentration of species A in the fluid 
(kmol m; °) 

dp particle diameter, equivalent diameter of sphere 
of the same external surface area (m) 

di internal tube diameter of the reactor (m,) 


Dae effective liquid diffusivity inside the catalyst at 
temperature T m? mp 17!) 

Dia liquid diffusivity inside the catalyst at temperature 
T (m mp E) 

Dao liquid diffusivity inside the catalyst at temperature 
To (m; m;! s7!) 


Eact activation energy for sodium borohydride hydrol- 
ysis on ruthenium catalyst (66,900 kJ kmol~ !) 

ftp two-phase factor (2.3 was used in current study) 

hfg heat of vaporization of water, assumed to be con- 


stant 2250 kJ kg7™! or 40,500 kJ kmol”! 

AHads heat of reaction for the adsorption of borohy- 
dride ion on the surface of ruthenium catalyst 
(—35,000 kJ kmol™!) 

AH;,xn heat of reaction for the sodium borohydride 
hydrolysis (—210,000 kJ kmol~!) 

kL reaction rate coefficient for 
Langmuire—Hinshelwood kinetic model 
(kmol kg cat! s7!) 


kı mass transfer coefficient from liquid to solid 
interface, based on concentration driving force 
(mẹ m;* s~!) 

kmo mass aaler coefficient for water vapor (m7!) 

K isotherm adsorption coefficient for borohydride 
ion on the surface of the catalyst (m? kmol~!) 

L length of the reactor (m) 

me mass flow rate of sodium borohydride solution 
(kg s~!) 

me initial mass flow rate of sodium borohydride solu- 


tion (kg s7 1) 
MWy, molecular weight of hydrogen (kg kmol” !) 
MWnņ,o molecular weight of water (kg kmol?) 


AB molar flow rate of liquid water (kmol s7 1) 

nc molar flow rate of hydrogen (kmol s~!) 

np molar flow rate of water vapor carried with hydro- 
gen stream (kmol s7!) 

Wiis total rate of water vaporization per unit catalyst 


mass (kmol water kg cat! s7!) 


esas q rate of water vaporization corresponding to the 
generation of hydrogen per unit catalyst mass 
(kmol water kg cat! s7!) 

Nevap.2 rate of water vaporization corresponding to 


mass transfer from catalyst surface to the 
bulk gas stream per unit catalyst mass (kmol 
water kg cat~! s7!) 


P back 


Psat,H,0(T) saturation pressure of water vapor at temper- 


P; 
AP fritz 
AProtal 
A 


Re, =dppius/u; Reynolds number of the liquid phase 


R 
RH 


Spellet 


Wiotal 


XA 
Xquality 


Z 


Greek letters 


E 
Es 


PH2 


PH2Ovapor density of the water vapor (kg m —3) 


D 


T 


Subscripts 

initial condition 
sodium borohydride 
bulk phase 

liquid water 
hydrogen 


backpressure of the reactor (bar) 


ature T (bar) 

total pressure of the reactor at location z (bar) 
total pressure drop across the fritzs (Pa, psi) 
total pressure drop across the reactor (Pa, psi) 
rate of reaction per unit catalyst mass 
(kmol kg cat™! s7!) 


universal gas constant, 8.314 (kJ kmol”! K7!) 
relative humidity of hydrogen stream 

the external surface area of the pellet (m?) 
temperature in the reactor at location z (K) 
ambient temperature (K) 

superficial velocity of the liquid phase through the 
bed (m s7!) 

superficial velocity of the gaseous phase (m s7!) 
average velocity of the multi-phase fluid across 
the reactor (m s7!) 

volumetric flow rate of the gaseous phase (m?) 
the volume of the pellet (m?) 

total mass flow rate of the liquid phase (kg s~!) 
total mass flow rate of the gas phase (kg s7!) 
total mass flow rate of the multi-phase fluid 
(kgs~') 

conversion of sodium borohydride 

fraction of the gas phase out of the multi-phase 
fluid 

axial direction (m) 


void fraction of packing in a packed bed (m my 3) 
internal void fraction of the catalyst pellet, 
between 0.3 and 0.8, typically 0.40 

thiele modulus number for cylinder 

overall effectiveness factor 

effectiveness factor 

liquid viscosity at temperature T (kg m7! s7!) 
liquid viscosity at temperature To (kg m7! s7!) 
catalyst bulk density in a packed bed (kg m~?) 
fluid density (kg m~?) 

density of the gaseous phase (kg m~?) 

density of catalyst (kg cat m; 3) 


density of the hydrogen gas (kg m~?) 


weighted average density of the multi-phase fluid 
(kgm~°) 
tortuosity factor, typically 3.0 to 4.0 
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Fig. 1. Section view of the reactor. 


stainless steel reactor. A distinct feature of this reactor is its 
stainless steel temperature profile probe, which enables accu- 
rate measurement of temperature profiles inside the reactor. The 
locations of all thermocouples (designated as PT-#) from the 
inlet are listed in Table 1. The inner diameter of the reactor was 
2.09 cm, and the height of the reactor was 28.7cm. The wall 
thickness of the reactor was 2.9 mm. The reactor had an internal 
volume of 98 mL. The catalyst used was 3% 2 mm Ru on carbon 
extrudate (Johnson Matthey). The total catalyst mass was 55.9 g, 
and the catalyst bed density was therefore 570 kg m°. At the 


Table 1 
Axial locations of thermocouples from the bottom of the reactor 


Name Distance from the 
inlet (cm) 
PT-1 28.70 
PT-2 25.53 
PT-3 22.35 
PT-4 19.18 
PT-5 16.00 
PT-6 12.83 
PT-7 9.65 
PT-8 6.48 
PT-9 3.30 
PT-10 0.13 


inlet, a stainless steel fritz was used to retain the catalyst and 
to distribute inlet flow. At the outlet, another stainless steel fritz 
was used to retain the catalyst. To simplify the physical model 
for reactor modeling, the reactor was thermally insulated dur- 
ing experiments to simulate adiabatic operation conditions. All 
downstream tubing and devices were also insulated for the same 
purpose. 

Tests were conducted at different fuel solution flow rates, 
operating pressures and inlet temperatures. The fuel concentra- 
tions used were 5%, 10%, 15% and 20%. The concentration 
of NaOH (stabilizer) in all solutions was 3% and the pH of 
all solutions was 14 [6,7]. The accuracy of temperature mea- 
surements was within +0.5 °C. Relative humidity was measured 
indirectly by condensing the hydrogen stream using a plate heat 
exchanger and measuring the weight of the condensate. The 
uncertainty of relative humidity is approximately +10% [3]. 
Conversion of NaBH4 was measured using two methods. The 
first involved measuring the amount of hydrogen gas evolved 
from acid-catalyzed hydrolysis of the product liquid. The second 
method involved measuring the flow rate of hydrogen generated 
using a flowmeter and comparing it to that based on the measured 
fuel inflow rate and 100% conversion. The two methods were 
found to match each other quite well, typically within 5%; as a 
result, only the second method was used for later experiments 
because it was easier to implement. 


3. Modeling 
3.1. Physical processes inside the reactor 


At the reactor’s entrance, liquid phase NaBHg solution (reac- 
tants) is introduced from the bottom of the catalyst bed. Then, 
the catalyst initiated hydrolysis that heats the reactants and the 
catalyst bed. The increase in reactor temperature along the flow 
path increases the reaction rate significantly. Because the reac- 
tion is exothermic (—210,000 kJ kmol7!), water vaporized along 
with hydrogen generation. Throughout the catalyst bed, the flow 
is multi-phase (liquid and gas) and multi-component (NaBH4, 
NaOH, NaBO2, H20, H2). Heat transfer involves convection 
and conduction, in addition to the heat conduction inside the 
catalyst pellets. All heat transfer processes are closely coupled 
with mass transfer processes (e.g., diffusion of NaBH; into the 
pores of the catalyst pellets) and chemical kinetics; and detailed 
modeling of all processes in the catalyst bed is therefore chal- 
lenging. This work reports a foundational effort in solving this 
complicated problem. 

The present model assumes that the catalyst is homoge- 
neously distributed inside the reactor, and that the reactor 
conditions vary in axial direction only and are uniform in each 
cross-section. For modeling purposes, the catalyst only affects 
the results through the void fraction of the catalyst bed. The 
kinetic model assumes uniform concentration and temperature 
at each cross-section and employs an overall effectiveness factor. 
The overall effectiveness factor accounts for both the external 
mass transfer coefficient outside the catalyst pellet and internal 
diffusion effects inside the catalyst pellet. A schematic of the 
model is shown in Fig. 2. 


J. Zhang et al. / Journal of Power Sources 170 (2007) 150-159 153 


NaBH,+NaBO,+NaOH+H,+ 
H,0(g)+H,O() 


N 


iii, Ue 


dz 


Insulation 


AAAA A A 


N 
N 
N 
N 
N 
N; 
X 
N 
Ñ; 
SS 
N 
N 
N 
N 
S 
SY 
N 
N; 
X% 
N 
N 
SS 
N 
N 
SS 
\ 
SS 
SS 
N 
N 
N 
N; 
X 
N 
N; 
X 
N 
N 
X 
N 
N 
X 
N 
N 
N 


YY Y@]@_ PY@E@EEE@q_uqux@@@u@PEaq@PE@@@uququqq@qB@@ 


Ve, 


z=0 
C, =C4) (NaBH, Concentration) 
T =T, (temperature) 
NaBH,+H,O()) ñp = Ny. (liquid water molar flow rate) 
+NaOH ; 


ñs =0 (H, molar flow rate) 
ñp = 0 (water vapor flow rate) 


Fig. 2. Schematic of the homogeneous model. 


3.2. Governing equations 


The concentration of component A (NaBHy) satisfies a 
steady-state continuity equation [8] (for simplicity, axial dif- 
fusion is ignored): 


Ms ie = rA Pb (2) 
where, 
me 
uls = (3) 
5 prAt 
pr = pro — xa) — prexa (4) 
Ca 
mal- A (5) 
Cao 
me = mfo — AC MWy, — tpMWy,0 (6) 


For other components, a molar flow rate approach is used. 
The steady-state continuity equations for component B (liq- 
uid water), component C (hydrogen) and component D (water 
vapor) can be written as: 


dig j j 

dz = (2r A Pb = NevapPb) At (7) 
die 

dz = —4r'y PoAt (8) 
din 

dz = NevapPoAt (9) 


P, Eam oT) Fa, 


Y Pasgl®) 


Fig. 3. Schematic of water vapor generation inside a catalyst pellet and water 
transport in the bulk space. 


where, 


/ 


Nevap = 


n 


evap, ir Nevap,2 (10) 


Eq. (8) for hydrogen is based on reaction stoichiometrics. The 
negative sign is present because the rate of reaction is defined 
based on the consumption of NaBH4. 

Modeling water vapor generation described in Eq. (10) is 
difficult because water vaporizes at different temperatures and 
locations in the reactor, in contrast to the isothermal conditions 
for most boiling situations [9]. Initially we assumed that water 
vapor generation is associated only with hydrogen generation 
and that hydrogen generated in each spanwise control volume 
was saturated with water vapor at the local temperature because 
hydrogen is generated from small pores inside the catalyst sup- 
port, as shown in Fig. 3. As a result, the partial pressure of 
the water vapor carried by the hydrogen would be the saturated 
water vapor pressure at that temperature. Therefore, according 
to Dalton’s law: 


—4r' Pb Z Pi — Psa,mo(T) 
Nevap, 1 Pb Prat, H,0(T) 


(1) 


The preliminary simulation results indicated that the pre- 
dicted relative humidity as defined below, of the exiting 
hydrogen stream was significantly lower than measured data. 

PH 0.g 
Psat, OCT) 


RH = (12) 
Therefore, the previous assumption was modified. Because of 
the large temperature variation inside the reactor, saturated 
hydrogen gas generated in the lower part of the reactor becomes 
unsaturated when it passes through the upper part of the reac- 
tor. Mass transfer of water from the liquid phase to the bulk 
gas stream should also be included. Mass transfer rates increase 
with an increase in driving potential (Psat — PH,0,g) and with 
the molar flow rate of the gas stream. As a result, we account 
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for additional mass transfer of water from the liquid phase to the 
bulk gas stream by relating it to the hydrogen molar flow rate. 


i kmo Psat — Puz0,¢ . 
Nevap,2 = NC 
PoAt Pt — Ph,o,g 


The denominator is used to convert the pressure-based driving 
potential to a molar flow rate-based driving potential using Dal- 
ton’s Law, similar to Eq. (11). The mass transfer coefficient for 
water kpo is not available in the literature and was estimated 
by varying its value to match one set of the experimental data 
as described in Section 4.1. In the calculations, P; is assumed 
to equal to the measured pressure head at the entrance of the 
reactor and remains constant throughout the reactor. 

Because the reactor was thermally insulated during exper- 
iments, the reactor is assumed to be adiabatic, and heat loss 
from the reactor to the ambient has been ignored. Therefore, the 
energy equation can be written as [8]: 


ůcMWn, Cp, Ho + ADMW4,0Cp,H0 dT 
At dz 


(13) 


(at, sens + 


= A Fyxnt', Pb a NégNevap Pb (14) 


To predict the pressure drop of multi-component, multi-phase 
flow through a porous medium, a homogeneous equilibrium 
model (HEM) from two-phase flow theory [10] may be used 
assuming that the liquid and gas phases have the same velocity. 
The total pressure gradient consists of three components (accel- 
erational, gravitational and frictional) that can be calculated as: 


dp _ dp dp dp 
dz total dz accelerational dz gravitational dz frictional 
(15) 

dp _ Wrotal da (16) 
dz accelerational Ay dz 

dp - 

- = jg (17) 
dz gravitational 
where, 
ae Wiotal (18) 

PAE 
Wiota = Wg + We (19) 
We = mf, — ANCMWy, — 7DMW3,0 (20) 
Wg = nAcMWy, + ApMWp, (21) 
W, 
Xquality = wiw (22) 
7 1 
p= (23) 
(Xquality / Pg) +- Xquality / Pf) 
CPH, + AD PHO va 
= 2 20 vapor 24) 


c + ñD 


For the frictional pressure gradient term in Eq. (15), no cor- 
relation is available for multi-phase flow through porous media; 
therefore, we propose the following formula based on single- 
phase flow through porous media [8] and two-phase flow theories 
[10]. 


dp 


_i 
, = 2 frep— (25) 


Z | frictional dp 

The equivalent diameter dp of the catalyst pellet (the diameter 
when the external surface area of the sphere equals that of non- 
spherical particles [11]) was used here instead of the reactor 
diameter d; to account for the porous-media flow and a two-phase 
factor frp was used to account for two-phase phenomena. We 
found that frp = 2.3 resulted in a good match with experimental 
data. 

To calculate the total pressure drop across the reactor, the 
pressure drop across both fritzes (A Pfritz) is determined by a cor- 
relation provided by the manufacturer. The total pressure drop 
can then be calculated as: 


L dP 
A Protal = aan 
0 dz 


The entrance (z=0) boundary conditions for the foregoing con- 
tinuity, energy and momentum equations are: 


dz + A Prritz (26) 


total 


Ca = Cao (27) 
T= (28) 
ñg = ÎB,0 (29) 
ic=0 (30) 
ip =0 (31) 


Chemical reactions are modeled using the Langmuir— 
Hinshelwood kinetic model developed in reference [2]. To 
account for the diffusion of reactants from the catalyst sur- 
face to the interior of the catalyst pellet, an effectiveness factor 
has been employed. The effectiveness factor calculation for a 
Langmuir—Hinshelwood kinetic model is quite complicated, and 
most prior results [12—14] have been presented on a log-log 
scale, making an integrated computer model difficult to imple- 
ment. As a result, we have used the general modulus and 
corresponding effectiveness factor expression proposed by Hong 
et al. [15]. The proposed modulus approximates the standard 
modulus to within 5%. To consider both the effects of mass 
transfer resistance from the bulk liquid phase to the surface of the 
catalyst pellet and internal diffusion effects, an overall effective- 
ness factor was used [11]. The primary equations for chemical 
kinetics are: 


KCa 


nen me 32 

TA NG LII KGa (32) 

ky, = AeW Pact/ RT (33) 
A Hads 

K = 1.6 x 1074 exp (-2) (34) 
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1 
= 35 
75 = On) + (kL Kps/kidg) aa 
“1 1 1 a 
"= % \tanhGd) 39 
_ Vpellet kL Kps (37) 
Spettet | Di A,e(2KC4 + (1/0 + KCa))) 


According to reference [2], the pre-exponential factor A in 
Eq. (33) is e!6903 kmol s~! kgcat., and the activation energy is 
66,900 kJ kmol~!. AHaas in Eq. (34) is —35,000kJ kmol7!, 

Two variables, the tortuosity factor t and internal void frac- 
tion €g, are needed to calculate the effective diffusivity (D),a ¢) of 
liquid inside the pores [8]. These factors account for the poros- 
ity of the catalyst support and the fact that the characteristic 
diffusion path length through the pores is greater than the mea- 
surable pellet thickness due to the tortuous nature of the pores 
and the pore constrictions. According to the recommendation of 
Satterfield [16], we have used a tortuosity factor of 4.0 and a 
void fraction of 0.4. The effective diffusion coefficient of liquid 
inside the pores is thus calculated from the intrinsic diffusion 
coefficient as: 


(38) 


To account for the significant variation of diffusivity with tem- 
perature, the following approximation has been adopted [11]: 


Mio T 
Di a = Diao 5 


m To a 


Dwidevi and Upadhyay [17] have developed a mass transfer 
coefficient model that is valid for both gases (Re > 10) and liq- 
uids (Re>0.01) in either fixed or fluidized beds. However, the 
mass transfer coefficient calculated using this model is of the 
order of 1076 m s—!, which is much smaller than mass transfer 
coefficients reported for similar conditions [11]. The correla- 
tions may not be applicable to multi-phase flow, which exists in 
our study. Because the generation rate is so much larger than the 
liquid flow rate, the actual mass transfer coefficient should be 
much larger. As a result, we have used a typical mass transfer 
coefficient value kj=2.0ms~! [11] and found that under this 
assumption, the overall effectiveness factor equals the internal 
effectiveness factor, indicating that the reaction is limited by 
internal diffusion. 

Eqs. (2)-(39) have been assembled and solved iteratively 
using Engineering-Equations Solver (EES) software [18]. The 
numerical integration grid spacing was 0.0005 m in all calcula- 
tions. Reducing the spacing further did not affect the results. 


4. Results and discussion 
4.1. Model parameters 

Both the diffusion coefficient Dj, 9 and the water vapor mass 
transfer coefficient ky,0 were adjusted to fit experimental data 


at a moderate flow rate (mf o = 48.5 g min” !) for 10% NaBH4. 
The density of the 10% NaBH; solution has been measured to 
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Fig. 4. Axial temperature and conversion profiles (fo = 48.5 gmin—!, 10% 
NaBHag, reference case). 


be 1045 kg m~?, and the density of the corresponding discharge 
solution has been measured to be 1200 kg m~*. The specific heat 
of the liquid phase is assumed to be that of water because of the 
relative low concentration of NaBHa. The heat of vaporization 
of water is 2250kJkg~!. The void fraction of the packed bed 
is 0.276, and the void fraction of the catalyst pellet is 0.4. The 
bulk density of the catalyst in the packed bed is 570 kg m~? 
while the density of the catalyst pellets is 640 kg m~+. The heat 
of reaction for NaBH; hydrolysis is —210,000 kJ kmol”! [19]. 
We found that Di ao =8.0 x 107!° m? s7! and kp,o = 3.0 m7! 
produced good fits. This value of Di A, is within the typical 
range of diffusion coefficients for liquids inside catalyst pores 
[20]. 

Figs. 4 and 5 present the predicted profiles for temperature, 
chemical conversion, relative humidity and molar flow rates of 
hydrogen and water vapor. The temperature increases slowly ini- 
tially. When it reaches 40 °C, the temperature begins to increase 
much more rapidly. Chemical conversion follows a similar trend 
and increases most rapidly when the temperature is between 
40 and 85°C. When chemical conversion approaches 100%, 
temperature peaks and then begins to decrease as a result of 
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Fig. 5. Axial relative humidity and molar flow rates profiles (mfo = 
48.5¢ min7!, 10% NaBHag, reference case). 
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evaporative cooling from the mass transfer of liquid water to 
the hydrogen stream. The profile of molar hydrogen flow rate 
follows the same trend as chemical conversion, as expected. The 
profile of molar water vapor flow rate generally follows the con- 
version trend except near the outlet, where the hydrogen flow 
rate remains constant and the water vapor flow rate continues 
to increase because of continuous evaporation of liquid water 
into the gas stream. Relative humidity initially decreases to a 
minimum and then begins to increase. The initial decrease in 
relative humidity is due to the fact that the water vapor added 
to the hydrogen is insufficient to match the increase in temper- 
ature because the temperature and reaction rate are low. The 
later increase in relative humidity is due to the opposite circum- 
stance, in which the added water vapor to the hydrogen stream, 
as a result of a much higher reaction rate, is more than enough 
to match the increase in temperature. 


4.2. Model predictions for 10% NaBH4 


With Dj,4,9 and kp,o obtained at one condition, temperature, 
chemical conversion, relative humidity and molar flow rates of 
hydrogen and water vapor profiles for other conditions were 
calculated and compared to experimental data [3]. 


4.2.1. Effects of flow rate 

Figs. 6—9 show the predictions together with experimental 
data under different flow rates. The predicted temperature pro- 
files (Figs. 6-8) match experimental data very well, as does 
the outlet chemical conversion except for the case of mo = 
64.0 g min—! where the predicted chemical conversion exceeds 
the measurement. The temperature profiles reveal that with 
increased flow rate, the predicted maximum temperature loca- 
tion (hot-spot) moves downstream, and the rate of conversion 
decreases because the reaction is highly sensitive to temper- 
ature. A relative humidity profile (Fig. 9) is only shown for 
mt = 64.0 g min™! because only a measurement under this con- 
dition was available. The relative humidity was underpredicted, 
and considering that the experimental uncertainty is approxi- 
mately +10%, the agreement in relative humidity is reasonable. 
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Fig. 7. Axial temperature and conversion profiles (nfo = 32.5 g min™!, 10% 
NaBH4). 
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Fig. 8. Axial temperature and conversion profiles (nfo = 64.0 g min—!, 10% 
NaBH4). 


4.2.2. Effects of inlet heating 

With inlet heating, temperature predictions near the inlet 
exceed measured temperatures slightly, but the relative humid- 
ity at the exit of the reactor matches experimental data very well 
as depicted in Figs. 10 and 11, respectively. As expected, the 
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Fig. 9. Axial relative humidity and molar flow rates profiles (m9 = 
64.0 g min™!, 10% NaBHy4). 
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Fig. 10. Axial temperature and conversion profiles (rng o = 64.0 g min—!, 10% 
NaBHy, inlet heating). 


relative humidity increased significantly from that in Fig. 9 at 
the same flow rate but without inlet heating. With inlet heating, 
the high-temperature portion of the reactor is longer, allowing 
more time for the transfer of water from the liquid phase to the 
gas phase while the gas stream passes through the reactor. 


4.2.3. Effects of operation pressure 

Under a higher system pressure condition (Fig. 12), the initial 
temperature rise is captured well by the model; however, the out- 
let temperature is overpredicted, as is the chemical conversion 
primarily due to the overprediction in temperature. Another pos- 
sible source of discrepancy in chemical conversion is the reduced 
contact between liquid fuel and catalyst at higher operation pres- 
sures as explained in Ref. [3]. This factor is not considered in 
the current model because the detailed mechanism for reduced 
contact between liquid fuel and catalyst is still unknown. 


4.3. Model predictions for 15% NaBH4 


Using the same parameters discussed in Section 4.1, the 
model was used to predict reactor behavior for 15% NaBH4 
solution at different flow rates. Because relative humidity was 
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Fig. 11. Axial relative humidity and molar flow rates profiles (9 = 
64.0 gmin—!, 10% NaBH, inlet heating). 
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Fig. 12. Axial temperature and conversion profiles (ngo = 64.0 g min—!, 10% 
NaBHa, 13.8 psig system pressure). 
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Fig. 13. Axial temperature and conversion profiles (go = 21.7 g min™!, 15% 
NaBH4). 


not measured in these tests, only temperature and chemical con- 
version profiles are presented. At lower flow rates, the model 
does not capture the initial temperature rise well, as shown in 
Fig. 13, while at higher flow rates the model captures the initial 
temperature rise very well, as shown in Fig. 14. As discussed in 
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Fig. 14. Axial temperature and conversion profiles (ngo = 64.3 g min™!, 15% 
NaBH4). 
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Table 2 
Comparison of predicted pressure drops with experimental data (10% NaBH4) 
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ring (g min—!) 


11.5 32.5 48.5 64.0 64.0, inlet heating 64.0, 9.6 psig 64.0, 13.8 psig 
Measured AP (psi) 0.8 2.0 3.0 3.2 4.3 2.4 1.8 
Predicted AP (psi) 0.90 22 3.0 3.2 6.3 29 2.8 


Ref. [3], the initial temperature rise was very different for flow 
rates below 35 g min™! than for flow rates above 35 g min™!. 
For flow rates below 35 g min™!, the rapid temperature rise may 
be partially due to axial heat conduction from downstream to 
upstream. The entry temperature increases to approximately 
40 °C (Fig. 13) while the fuel’s initial temperature is only 22°C. 
With an increased fuel flow rate, the effects of heat conduction 
are reduced, and the model predictions agree with experimental 
data. Temperatures near the exit of the reactor are underpredicted 
for all cases, presumably because of overestimated water evap- 
oration into the hydrogen stream passing through the reactor. A 
smaller value of ky,0 may mitigate this situation. Further, using 
a smaller ky,o may be consistent with the experimental process 
because at higher NaBH4 or NaBO? concentrations, water is 
less likely to evaporate than it is at lower concentrations [21]. 
The chemical conversion predictions match experimental data 
well for flow rates up to 53.2 g min~!. At the higher flow rate of 
64.3 g min™!, chemical conversion is slightly overpredicted. 


4.4. Prediction of pressure drop 


The predicted pressure drops across the reactor are com- 
pared to experimental data in Table 2. Only experiments with 
10% NaBH; are considered. At 15% NaBHy concentration, 
the viscosity of the discharge stream increases significantly 
as explained in reference [3], and the pressure characteristics 
change dramatically. The pressure drop predictions are excel- 
lent, especially for different flow rates. The predicted pressure 
drop under the inlet heating condition exceeds the measurement, 
perhaps due to the overpredicted initial temperature rise dis- 
cussed in Section 4.2.2. For pressurized conditions, the pressure 
drop is also overestimated, perhaps due to the discrepancies in 
chemical conversion. However, the trend of decreasing pressure 
drop at higher system pressures was captured by the present 
model. Decreased pressure drop is likely due to lower water 
vaporization rates inside the reactor at higher operating pres- 
sures. 


4.5. Recommendations 


Both the previous experimental study [3] and present mod- 
eling results reveal that inlet heating can significantly increase 
reaction rates in the entry section of the reactor and increase 
the overall chemical conversion at high flow rates. Inlet heat- 
ing (e.g., via exhaust recirculation) could be adopted to reduce 
the length of the reactor and the amount of catalyst used. Fur- 
ther, temperature measurements could prove to be effective in 
real-time performance (conversion) monitoring. The effects of 


operation pressure on overall chemical conversion were incon- 
sistent between experimental study and modeling results, and 
further investigation of this issue is warranted. 


5. Conclusions 


A one-dimensional/homogeneous catalyst distribution model 
has been developed to simulate thermo-chemical processes 
inside a NaBH; hydrolysis packed bed reactor. The model suc- 
cessfully captures water vapor generation both by chemical 
reactions and by mass transfer from the liquid phase to the gas 
phase. The pressure drop has been calculated based on two-phase 
flow and porous media theories. The overall model has been used 
to predict reactor behavior at different fuel concentrations, flow 
rates and reactor pressures. The predicted temperature profiles, 
fuel conversion and relative humidity compare reasonably well 
to experimental data. The predicted pressure drops also match 
experimental data very well at different flow rates, although pres- 
sure drops under pressurized and inlet heating conditions are 
overpredicted. Considering the complexity of processes occur- 
ring inside the reactor, the present results are encouraging. The 
numerical model is based on fundamental heat and mass transfer 
and chemical kinetics theories; therefore it is expected to be use- 
ful in studying and optimizing hydrogen storage systems using 
NaBH; and other chemical hydrides. 
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